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Abstract. We consider the random continued fraction 



S(t):= , te 



t 

S2 



S3 + • ■ • 

where the independent random variables with the same gamma distri- 

bution. Every realisation of the sequence defines a Stieltjes function that can 
be expressed as 

io l + xt 

for some measure a on the positive half-line. We study the convergence of 
the finite truncations of the continued fraction or, equivalently, of the diagonal 
Fade approximants of the function S. By using the Dyson— Schmidt method for 
an equivalent one-dimensional disordered system, and the results of Marklof et 
al. (2005), we obtain explicit formulae (in terms of modified Bessel functions) 
for the almost-sure rate of convergence of these approximants, and for the 
almost-sure distribution of their poles. 



1. Introduction 

Let s — (si,S2,...) be a sequence of positive real numbers and consider the 
analytic continued fraction 

(1.1) S{t):^ ^-^ , teC\M_. 



Sl 



t 

S2 



S3 + • • • 

This continued fraction defines a Stieltjes function; it can be represented in integral 
form as 

for some measure a supported on the non-negative half-line such that the moments 

POO 

(1.3) m„ := / x"cr(da;), n e N, 

Jo 
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exist. By an obvious use of the geometric series, every Stieltjes function can be 
expanded formally in powers of t: 

oo 

(1.4) S{t) - m-j {~ty as t 0+ . 

3=0 

Hence 5* is the moment generating function of the measure a. It is a well-known 
fact of great practical importance that, given the first n of the moments, one may 
construct the rational function 

(1.5) ' ^^^\«-' 

si H 



S2 



t 



where 



I n/2 — 1 if n is even , , ^ I '^/2 if n is even 

degP„ — < and desQr, = < 

[(n-l)/2 if n is odd \{n-l)/2 if n is odd 

This truncation of the continued fraction has a MacLaurin expansion whose 
nth partial sum agrees with that of the series ()1.4p . Hence Sn is a diagonal (if n is 
odd) , or near-diagonal (if n is even) Pade approximant of S. 

Now suppose that the s„ are independent positive random variables with the 
same distribution, say /z. We shall consider the following questions: 

(1) What are the almost-sure analytic properties of these Stieltjes functions? 

(2) What is the almost-sure leading asymptotic behaviour of the error S{t) — 
Sn{t) as n ^ oo? 

These questions are of interest because Pade approximation is widely used in 
applied mathematics as a practical means of accelerating the convergence of the 
partial sums of series obtained by perturbation methods. As pointed out by Ben- 
der & Orszag (1978), the consideration of many particular cases where the s„ are 
deterministic reveals a wide range of large-n behaviours. Our motivation for study- 
ing the random case is to gain some insight into the asymptotic behaviour of Pade 
approximation in the "generic" case. 

Our study of diagonal Pade approximation reduces to aspects of the large-n 
behaviour of the denominators Q„. For this reason, as in the deterministic case, 
the cornerstone of the analysis is the three-term recurrence relation 

(1-6) Qn+l = tQn-l + Sn+lQn , Q-1 — 0, Qo — ^ ■ 

(The Pn satisfy the same recurrence relation, albeit with different initial condi- 
tions.) This recurrence relation makes a link between Pade approximation and a 
rich set of other mathematical entities, such as orthogonal polynomials, products 
of random matrices, and discrete Schrodinger-like operators. By exploiting results 
that are well-known in these related fields, one may obtain — for a very large class of 
distributions /i of the coefficients s„ — some partial answers to the questions stated 
earlier. Our contribution in the present paper is to elaborate the particular case 
where /i is the gamma distribution. More precisely, we obtain explicit formulae (in 
terms of Bessel functions) for the leading term in the asymptotic behaviour of the 
error of Pade approximation and for the asymptotic distribution of the poles, as 
well as the location of the essential spectrum of the measure a. 
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In the remainder of this introductory section, we describe briefly the key ideas 
underlying the analysis. Then we summarise our main results in the form of a 
theorem. 

1.1. The moment problem. The Stieltjes moment problem is, given a sequence 
{'T^njneN, to determine whether or not there exists a measure a such that equation 
(jl.3p holds for every n. Historically, mathematical objects such as the analytic 
continued fraction (jl.ip . orthogonal polynomials and Pade approximants were in- 
troduced as tools in the study of this moment problem (Akhiezer 1961; Nikishin & 
Sorokin 1988; Simon 1998). Stieltjes (1894) showed that a necessary and sufficient 
condition for the existence of a measure a with the prescribed moments is 

(1.7) VneZ+, s„>0. 
He also showed that 

oo 

(1.8) s« = oo 

n=l 

is a necessary and sufficient condition for the uniqueness of the measure. 

Let us assume that condition (|1.7|1 holds and describe in very broad terms one 
way of "reconstructing" a from its moments; see Akhiezer (1961) and Nikishin & 
Sorokin (1988) for a detailed treatment. 

Recall that x' is a point of increase of the measure a if 

Ve>0, / cr(da;)>0. 

max{x' — 6,0} 

The spectrum of cr is the set of its points of increase and will be denoted spec(cr). 
For a; > 0, we shall denote by 5x the probability measure on R+ whose only point 
of increase is x. 

Given the m„, we may define an inner product, say (•, •)ni, on the space of 
polynomials as follows: if p and q are two polynomials with coefficients pi and qi 
respectively, then 

{P,q)ra := ^m^+jPiqj . 

Knowing the moments, we may compute the s„ and the S„. We remark that P2n 
and Q2n are polynomials of degree n — 1 and n respectively. Set 

(1.9) ^„(A) := ^^^\"Q2n{-l/\) . 

Then the fact that S2n matches the moment generating series to 0{t^^) implies 
that ipn is orthogonal to every polynomial of degree less than n, in the sense of the 
inner product (•, ■)m- It follows (see Akhiezer 1961, Chapter 1) that the roots of 
ipn are simple and lie in R+; denote them by 

< A„4 < Xn,2 < ■ ■ ■ < Xn,n < OO . 

Gaussian quadrature then defines a discrete measure 

n / n— 1 

(1.10) (T„ ;= ^cr„j(5A„.j , where (t„j = ^ Vl(^«j) 

3=1 \e=o 
that converges weakly to a measure a that solves the moment problem. In partic- 
ular, the inner product (•, •)in coincides with the inner product in L^(R+). 
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The moment problem can also be approached from the point of view of operator 
theory. The recurrence relation (ll.6|l for the Qn implies the following recurrence 
relation for the ^n'- 



where 



/l„_lV'n-l + Vn^n + ^nV'n+l 



if n 
if n G 




A-^o if n = , 
AiAn if n G Z_|_ 

1 



S2n+2\/S2n+lS2n+3 



n G N. 



These numbers may be used to define a certain Jacobi operator, say ^ , with a 
domain contained in the Hilbert space ^^(N), as follows: first, we consider sequences 
^ = (^0, Cij C21 • ■ •) with only finitely many terms and set 



(1.11) 



Wo Co + /loCi 

hn-l£,n-l + W„C„ 



if n = 
hn^n+i otherwise 



Given condition (|1.8p , it is then possible to extend this definition uniquely to obtain 
an essentially self-adjoint operator; we use the same symbol ^ to refer to this 
extension. It may then be proved that the moments of the spectral measure of the 
operator ^ are precisely the m„, and so this spectral measure coincides with tr; 
see Nikishin & Sorokin (1988). 

1.2. The density of states. Consider the finite-dimensional truncation 
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of the operator ^ . The spectrum of ^„ is the set of zeroes of the polynomial 
ipn defined by equation (|1.9p . By the Chebyshev-Markov-Stieltjes theorem (see 
Nikishin & Sorokin 1988, §2.8), between any two zeroes of there is a point of 
increase of a; so we are led to studying the distribution of the A„j . 
Define a measure k„ on R_|_ by 



(1.12) 



1 



This measure is the normalised eigenvalue counting measure of the matrix 
Indeed, we have 



iV„(A) 



# {j ■■ K,j < A} 



K„(dA') 



JO 

The normalised counting measure Kn has a weak limit, say k, as n — > cx), and so 
there is a function iV, called the integrated density of states of ^ , defined by 

^(A) / K(dA) = lim 7V„(A) . 
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If K is absolutely continuous, one can also speak of the density of states, say g, 
defined by 

(1.13) K(dA) = e(A)dA. 

Although the measures k and a may be very different, their essential spectra are 
the same. In the context of Pade approximation, the integrated density of states 
describes the distribution of the poles of the approximants. 

1.3. Krein's string. There is an interpretation, due to Krein, of the spectrum 
of the operator ^ in terms of the characteristic frequencies of a vibrating string 
(Akhiezer 1961). Consider a weightless, infinite, perfectly elastic string, tied at one 
endpoint x = 0, along which some beads are distributed. Let S2n be the mass of 
the nth bead, and denote by (xn,yn) its position in the xy-plane. We assume that 
the fixed and given by the recurrence relation 



S2n+1 , = . 



For a string of uniform unit tension, the small vertical motion is then described by 
the discrete wave equation 

(1.14) S2nyn = , n€Z+ . 

S2n+1 S2n-1 

To study the characteristic frequencies of the string, we set 

Vn+l - Vn 



Vn = Vn sin {tot) and £,„ 

S2n+1 



Then equation (|1.14p reduces to 

1 



S1S2 S1S2 

and 



6 Co = -t^^-^o 



"Cn+l ( ^ I ^ S,n-1 ~ — W^^„ 



S2n+lS2n+2 S2n+1 \S2n S2n+2 J S2riS2' 

where n> 1. Comparing this with the definitions of ^ and given earlier, it is 
readily seen that 

, (-l)"go , , 2^ 
U = iPn{UJ ) . 

VS2n+l 

1.4. The complex Lyapunov exponent. Dyson (1953) developed a method for 
studying the characteristic frequencies of the one-dimensional disordered chain 

(1.15) j/n = C2„-l(j/„+l - y„) - C2„-2(2/« - J/n-l) , n £ Z+ . 

Here, C2n-i and C2n-2 are the ratios of the elastic modulus of the nth. spring and 
of the mass of the two particles attached to it. Disorder may be modelled in 
many ways; for instance by assuming that the c„ are independent and identically 
distributed. The approach was later simplified by Schmidt (1957) and applied to 
the tight-binding Anderson model for a one-dimensional crystal with impurities. 

Luck (1992) gives a very readable, well-motivated account of the Dyson-Schmidt 
approach; in brief, it builds on the intimate connection between second-order dif- 
ference equations, continued fractions and Markov chains. For our purpose, it will 
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be convenient to work with the random difference equation 

(1.16) u„+i - M„_i = ^^u„ , n = 0,l,2,... 

where t is a parameter in C\M_ , and \f ■ is the branch of the square root function 
defined on C\]R_ that returns a number with a non-negative real part. We note 
the obvious 

Lemma 1.1. For every t G C\M_ . 

Qn{t) = (^Vi^ Un , 

where Un solves the difference equation I11.16\) with it_i = and uq = 1. 
The relevant continued fraction is 

(1.17) Z:^ViS{t)^ 



1 



1 



S2 



and we write 

(1.18) Zn^VtS^ 

for its truncation. Let u-i and uq be complex random variables. Then equation 
(|1.16p defines a sequence of general term u„ by recurrence. The distribution j/^ of 
the random variable Z is a stationary distribution for the Markov chain 

(1.19) Z = (Zo, Zi, Z2, ...) , where Z„ := ^ . 

(In the terminology of iterated random maps, the random variables Z„ and Z^ 
are, respectively, the backward and forward iterates associated with the continued 
fraction; when m_i = and uq = 1, they have the same distribution, but their 
asymptotic behaviours are very different; see Diaconis & Freedman 1999.) 

It follows that the growth of the u„ may be quantified by means of the complex 
Lyapunov exponent defined by 

(1.20) A^{t) :=- [ Inziy^(dz), 

Jc 

where In . denotes the principal branch of the logarithm. Indeed, standard results 
from the theory of Markov chains imply that, if Z has a unique stationary distri- 
bution, then 

In u„ In Un 1 x-^ , A . / N 

(1.21) ^ = ^--^InZ, >A^{t) 

for almost every realisation of s, independently of the choice of Zq (Breiman 1960; 
Furstenberg 1963; Meyn & Tweedie 1993). In particular, we have the formula 

lini ^ Rc [AM)] . 

n— >cx> Tl 
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Equation ()1.2ip is also central to the study of the integrated density of states of the 
operator ^ introduced earlier (Dyson 1953; Schmidt 1957; Luck 1992). We shall 
see that, under a very mild assumption on the distribution /i of the s„, 

7V(A) = --Im[A^(-l/A + iO+)] . 

TT 

1.5. Furstenberg's theorem. In order to carry out this programme, we shall also 
make use of the connection between the Markov chain Z and the product of random 
matrices 

(1.22) Un:=A^An-x---Ax, n = 0,l,2,..., 

where 



(1.23) A. 



1 

1 ^ 



The distribution /i from which the s„ are drawn induces, via equation p.23p . a 
distribution /2 on the group of unimodular 2x2 matrices. The fundamental results 
of Furstenberg & Kesten (1960) and Furstenberg (1963) — which are commonly 
referred to as 'Furstenberg's theorem' — imply in particular that, under very mild 
assumptions on /i, there is a unique measure vp_ on the group of unimodular matrices 
that is invariant under the action of the random matrix (jl.23|) . Furthermore, the 
number 

(1.24) 7^ ^E(lnK,|), 

which quantifies the growth of the product Un and is independent of the choice of 
matrix norm | • |, may be shown to be strictly positive. These results are of great 
relevance to our problem for two reasons: firstly, any invariant measure f/^ yields a 
measure Vfj^ that is stationary for the Markov chain Z, and vice- versa; secondly, 

7^ = Re [A^(t)] . 

Hence we deduce at once the uniqueness of the measure , as well as the exponen- 
tial growth of the w„: 

Proposition 1.1. Let the Sn he independent random variables in with a com- 
mon distribution /i that has at least two points of increase. Let {u„}ngN be the 
sequence defined by the recurrence 11.16\) . Suppose that 

s'^ /i(ds) < oo 

for some e > 0. Then, for almost every realisation of the sequence s, the following 
holds independently of the starting value Uq ^ 0; For Lebesgue-almost every t G 

lim = AM) and Re [AJt)] > . 

n — *oo 77, 

Furthermore, if we set t = —x + iO±, then for Lebesgue-almost every x > 0, 
In u 

lim -^AJ-x + iO±) and Re [A„(-a; + iO±)] > . 



The proof of this proposition is provided in Appendix O 
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1.6. The gamma distribution: statement of the main result. By using such 
machinery, we are able, for a very wide class of distributions fi, to deduce the almost- 
sure exponential nature of the convergence of diagonal Pade approximation and also 
to deduce the almost-sure singularity of the measure a. A more quantitative study 
requires the calculation of the complex Lyapunov exponent, but there are very few 
known instances where it can be expressed in terms of familiar functions. 

In his seminal paper on the disordered chain (jl.isp , Dyson studied in some detail 
the particular case where the c„ are independent and gamma-distributed. Dyson 
found the invariant distribution of the continued fraction 

Cot 



Cit 



1 



C2< 



1 



in the particular case where t > 0; he then used analytic continuation to obtain 
an expansion for the complex Lyapunov exponent at f < 0, and hence for the 
distribution of the characteristic frequencies. Dyson's continued fraction is not 
equivalent to ours (compare equations (jl.lSp and (|1.14p ). and so his analytical 
results do not transfer to our problem. However, the continued fraction p.l7p with 
independent, gamma-distributed s„, i.e. where 

(1.25) n{Ax) ^ x°-~'^ exp(-g;/6) dx , a, 6 > . 

was examined also by Letac and Seshadri (1983); they obtained the probability 
distribution and found an explicit formula for the corresponding Lyapunov ex- 
ponent for < > 0. In a recent paper, we generalised this result by finding and the 
real part of the complex Lyapunov exponent for every complex t (see Marklof et 
al. 2005); a straightforward extension of these calculations leads to the remarkably 
simple formula 



(1.26) 



Ka 



In this expression, da denotes differentiation with respect to a, and Ka is the modi- 
fied Bessel function of the second kind (see Appendix|B]). The following summarises 
the key results of this paper. 

Theorem 1. Suppose that the s„ are independent draws from the gamma distri- 
bution with parameters a > and b > 0. Then, for almost every realisation of the 
sequence s, the following holds: 

(1) The density of states is given explicitly by the formula 



^^^^ = 21^'^ 



1 



(2) spec(cr) — [0, oo) and its absolutely continuous part is empty. 

(3) For Lebesgue- almost every t G 



lim 

n — >oo 



\n\S{t) - Sn{t)\ 



= -2a„ln 
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1.7. Relation to other work. As stated earlier, our focus in the present pa- 
per is on the performance of diagonal Pade approximation, viewed as a method 
of summing some random series. Related questions have been considered in the 
past, in different contexts. Foster & Pitcher (1974) study the convergence of ran- 
dom T- fractions; these are continued fraction expansions which are in a one-to-one 
correspondence with the space of formal power series, but whose convergents are 
not Pade approximants. Foster & Pitcher (1974) show that, under very general 
conditions on the distribution of the coefficients, the difference between two suc- 
cessive convergents tends to zero exponentially fast, and that the exponent is twice 
the Lyapunov exponent associated with an infinite product of random matrices. 
Geronimo (1993) studies the random measure (on the unit circle) generated by a 
three-term recurrence relation with random, identically distributed coefficients; he 
shows the positivity of the corresponding Lyapunov exponent and deduces that the 
random measure is singular with respect to the Lebesgue measure. This list is not 
exhaustive; see also Csordas et al. (1973) and Mannion (1993). 

The question of the nature of the measure a has a counterpart in the theory of 
disordered systems which has been studied extensively in the context of Anderson 
localisation. For example, the tight-binding Anderson model uses a discretised 
version of the Schrodinger equation with a potential that takes random, identically- 
distributed values at every point in a doubly-infinite lattice. The resulting operator 
has a second-order finite-difference form like that of the operator ^ — in which, 
more precisely, /i„ = 1 and the w„ are independent and identically distributed — 
but it acts on sequences in (Z) . For a very wide choice of the distribution of the 
potential values w„, the Lyapunov exponent of the discretised Schrodinger operator 
is strictly positive, so that, by Ishii's formula, the absolutely continuous spectrum is 
empty. A more refined study (see for instance Carmona & Lacroix 1990 and Pastur 
& Figotin 1992) reveals that these operators have a pure point spectrum — that 
is, the spectrum is the closure of the discrete spectrum. Such operators are said to 
exhibit the localisation property because, for this type of spectrum, the generalised 
eigenfunctions decay exponentially fast as |n| — > oo. The rigorous extension of 
such detailed results to the semi-infinite case would involve technicalities which are 
outside the scope of the present paper. 

Our analysis exploits a number of ideas and techniques found in these earlier 
studies. We view our main contribution as that of exhibiting an interesting example 
of a class of random Stieltjes functions for which the leading behaviour of the error 
of diagonal Pade approximation, and the density of states of the corresponding 
Jacobi operator, are given explicitly in terms of special functions. 

The remainder of the paper is devoted to a detailed proof of theorem [TJ the 
first statement follows immediately from Dyson's formula for the density of states, 
which is derived in f|2] In f|3l we deduce the singularity of the measure a from 
the positivity of the real part of the Lyapunov exponent. In 21 show that the 
error of diagonal Pade approximation is inversely proportional to the square of the 
Un] this yields the third statement in the theorem. Finally, in ^ we provide a 
numerical illustration of our results. 

2. The formula for the density of states 

So that we can use proposition 11.11 we shall henceforth suppose that the s„ are 
independent draws from a distribution /i on R_|_ such that 
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(1) /i has at least two points of increase. 

(2) There exists £ > such that 




s^^(ds) < oo . 



To avoid needless repetitions, we shall not mention these particular assumptions 
explicitly again in the statement of the intermediate results leading to theorem [2l 
We begin by relating the growth of the 4'n to the complex Lyapunov exponent. 

Lemma 2.1. For almost every realisation of the sequence s, we have, for Lebesgue- 
almost every X G C\M+, 

li,^ Wa^=i, + 2A,(-1/A) 
and, for Lehesgue- almost every X G 

hm ^^^r.{X + iO±) ^.^^ ^ .^^^ 

n — > OQ 77, 



Proof Let A G C\K+. By definition, 
Hence, by lemma [TTTl 

/ s 2n 

IpnW = V^i^A" (y^-l/Xj U2n = (-l)"VS2n+l "2n , 

where m„ solves the difference equation (II. 16^ with t = — 1/A. This yields 
ln'0„(A) lns2„+i . , „lnM2n 

= h ITT + 2 . 

n 2n 2n 

The first statement in the proposition then follows from corollary lA. 21 The proof 
of the second statement is identical. □ 

Next, we examine the implications of the lemma for the distribution of the A„.j. 
By virtue of the recurrence relation satisfied by the ip„, we can write 

/ 2" \ 

■!/'n(A) = V*2n+1 I (A - A„,i) • • • (A - A„_„) . 

\j = l J 

Let 

En :== {A„j ■.l<j<n} 

and let X ^ En- Then 

lni/;„(A) ln.S2„+i 1 



fry-,^ Ini/'ri A ln.S2„+l 1-^ 1-^ 
2.1 = — + - > Insj + - > ln(A- A„j). 

j=i i=i 

Proposition 2.1. Suppose that 

/ |ln s| /i(ds) < oo . 
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Then, for almost every realisation of the sequence s, for Lebesgue- almost every 
X e C\R+, 

ln|A - A'l K«(dA') > 2Re[A^ (-1/A)] - 2 / lns^i{ds). 







Proof. Take the real part in equation (|2.ip . Then 

ln|i^„(A)| 



poo 1 

J ln|A-A'|Kn(dA') = -5]ln|A-A„,,| = 



1 1 2n 

lnS2„+l 1 ■r;^ 

— > ins," 



By lemma [131 the first term on the right tends to 

2Rc[A^ (-1/A)] , 

the second term tends to zero and, by the ergodic theorem, the third term tends to 

2 / Ins^(ds). 

□ 

Corollary 2.1. Under the same assumption, for almost every realisation of the 
sequence s, the sequence {Kn\n&i has a weak limit, say n, which is a probability 
measure on R+ . In particular, 

lim iV„(A) =N{\) := ( k{A\') . 







Proof. The proof is a specialisation of that given by Goldsheid & Khoruzhenko 
(2005) in the more general case of a non-Hermitian Jacobi matrix. See Appendix 
O for the details. □ 

Theorem 2. Let the s„ be independent random variables in R-|_ with a common 
distribution fi that has at least two points of increase. Suppose also that 

/ |lns| n{ds) < OO and / /J,(ds) < oo 

for some e > 0. Then, for almost every realisation of the sequence s, for Lehesgue- 
almost every A € M+j 

7V(A) = --Im[A^(-l/A + iO+)] . 
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Proof. Let A e M+\£;„. Then 



2n 



In ^„ (A + iO±) - 1^1^!^ " ^ H 1^ = 51 - ^" J + iO±) 

= ^ ln(A - A„j + iO±) + ln(A - A„.j + iO±) 

A„.j<A A„,j>A 

= ^ ln|A-A„j|+ (±i7r + ln|A- A„j|) 

A„j<A A„j>A 



^ln|A-A„j|± ^ iTT. 

J = l A„,j>A 



Hence, we have the identity 



lni/'n(A + iO±) _ hiS2K+i _^ }_ 
n 2n n 



2n 



El- 



(2.2) 



Take the imaginary part. The rcsuh then follows from lemma 12.11 and corollary 
O □ 



poo 

+ / ln|A-AVn(dA')±i^[l-iV„(A)] 
Jo 



Corollary 2.2. Suppose that the Sn are independent and gamma-distributed with 
parameters a and b. Then, for X > 0, we have the following formula for the density 
of states: 



TT^A 



1 



Proof. Let A > and set 



b^/x' 



For t = — 1/A + iO+, we have 



Hence 



and so 



\ b J 



daKa 



l"r]=-^2^4"r]-2' 



[daYa{z)+idaJa{z)] 



^a(¥ 
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We deduce the formulae 



(2.3) 
and 
(2.4) 



Im 



Re 



A„ 



iO- 



iO- 



^ Jaiz)daJa{z)+Ya{z)daYa{z) 
Ja{zY+Ya{zr 

TT , Ya{z)daJa{z) - Ja{z)daYa{z) 



2 Ja{zY+Ya{zf 

The proposition, together with the last equation, then yield 



1 - iV(A) 



2Ya{z)daJa{z) ~ Ja{z)daYa{z) 
Ja{zY+Ya{zf 



Differentiating both sides with respect to A, we find 



2 dz d 
TT dA dz 



da 



arctan 



Ya{z) 



Ja{z) 



1 



TT 6A3/2 dz 



'-da 



da 



arctan 



arctan 



Ya{z) 



Ja{z) 



We obtain the desired result by changing the order of differentiation on the right- 
hand side, and making use of the identity 

2 

■KZ 

□ 



Ja{z)Y:,{z)-Ya{z)j'M 



Corollary 2.3. Under the same assumption, for almost every realisation of the 
sequence s, 

spec(cr) = [0, oo) . 



Proof. By differentiating the identity (see Watson 1966, §13.73) 

8 

Jliz) + Y^{z) = ^ / Kn (2zsinht) cosh(2at) dt 
I" Jo 

with respect to a, we deduce that g is strictly positive. □ 

3. Singularity of the spectrum 

Corollary 12 .31 implies in particular that the radius of convergence of the generat- 
ing series of the moments is zero almost surely; in other words, the random Stieltjes 
functions that we have constructed are not analytic at the origin. 

The problem of determining the nature of the spectrum is more delicate. Every 
measure may be decomposed into three disjoint parts: its absolutely continuous, 
singular continuous and discrete parts, denoted by aac, o^sc and ad respectively. Ishii 
(1973) and Yoshioka (1973) showed that the spectrum of the absolutely continuous 
part is given by the formula 



(3.1) 



spec (ctqc) = < a > 



lim RcA^(t)=0 



This result may be established by examining the resolvent of the Jacobi operator 
^ . The following result is essentially equivalent. 
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Proposition 3.1. Let the Sn be independent random variables in M+ with a com- 
mon distribution fi that has at least two points of increase. Assume that 

/ |lns| /^(ds) < oo and / s^/i(ds) < oo 

Jr+ Jk+ 

for some e > 0. Then, for almost every realisation of the sequence s, 

spec(CTac) = ■ 



Proof. By lemma 12.11 and corollary IA.4[ for almost every realisation of s, for 
Lebesgue-almost every A G K+, we have 

(3.2) limi!i^>0. 

n — ^CJO ri 

Now, let > and consider the set 

{A e M+ : ?A„(A) = a {y^[\nn]^+'^) as n ^ oo} . 

By equation p.2p . this set has Lebesgue measure zero almost surely. On the other 
hand, it follows from the Men'shov-Rademacher theorem (see Nikishin & Sorokin 
1988, proposition 8.3) that for cr-almost A e R+, 

^/;„(A) = o (^/n [Inn]^^'') as n — > cxj . 

So, almost surely, for every cr-measurable set A, 

a(dA)= /" <j{dX)= f aacidX)+ f (T,,(dA) + /" ^^(dA) 
A JAns,, JAnSjj JAns„ JAnSr, 



asc{dX)+ / fT<i(dA)=: / (Tsc{d\)+ / crd(dA) 
Ans„ JAnSr, J A J A 



4. The rate of convergence 



□ 



Let be the sequence defined by the recurrence relation p.l6p with the 

starting values u_i = and uo = 1- Also, denote by £/' the shift operator on the 
space of complex sequences, i.e. 

3^ (Si, S2, • ■ •) = (S2, S3, . . .) . 

In order to emphasise the dependence of the continued fraction (|1.17p on the se- 
quence s, we shall sometimes write Z{s) and Z„(s) instead of Z and Z„. We have 
the following convenient representation of the error: 

Lemma 4.1. 



('_1^r^ / " 

z(s)-z„(s) = ^^ \{z{,r^^ 



Proof. By using the recurrence relations (|1.6[) and the identity 

. _„ , 1 
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it is straightforward to show (by induction on ri) that 

1 QnS Pn 



Z(^"s) = -- 



\/i Qn-lS — Pn-1 

Then 



3=0 



Lemma [TTT] then yields the desired result. □ 

Theorem 3. Let the Sn be positive independent random variables with a common 
distribution that has at least two points of increase. Suppose also that there exists 
e > such that 



^{As) < oo . 

^+ 

Then, for almost every realisation of the sequence s, for Lebesgue- almost every 
t G C\R_, 

limi^i^W^M^.-2Re[A,W]. 

n — >oo Ji 



Proof. Let t e C\M_ be fixed. We have 

S{t) ^ VtZ{s) and Sn{t) = VtZ„{s) . 
Hence, by lemma 14711 



\n\S{t)-S.^it)\ \n\t\ InKI ^ ly^i 

Consider the almost-sure limit of each of the three terms on the right-hand side 
of this equality as n — > cx3. The first term tends to zero. By proposition II. li the 
second term tends to 

-Re[A^(i)] . 

Finally, it follows easily from the ergodic theorem that the third term tends to 
the same limit. The result follows from a standard argument involving the use of 
Fubini's theorem. 

□ 

5. A NUMERICAL ILLUSTRATION 

Following the example of Dyson (1953), it is instructive to begin with an exam- 
ination of the (determistic) case where 

Set 

S°°{t):= i , ieC\M_. 

1 + 

t 
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Then 



1 + VI + 4t Jo 1 + xt 
where 

j^\/4:/x- Idx if < a; < 4 



a°°(dx) 



if X > 4 



Note that fi°° has only one point of increase, and so the hypothesis of Proposition 
13.11 is not satisfied. Indeed, for this choice of fi, the spectrum of the measure a is 
absolutely continuous. 

In this case, the continued fraction (|1.17p reduces to 



^l/^ + 4+l/^/^ 
and so the complex Lyapunov exponent is 



A^oo(t) = - / In z i^^^. (dz) = - / In z (5zoc (dz) = In — ^ ^ ^ 



In particular, an elementary calculation shows that 



A^°°(A) = 



1 - - arccos ^ if < A < 4 
1 if A > 4 



Next, let a G Z+ and denote by /i" the gamma distribution with = 1/a. As 
Dyson remarked, this distribution has mean 1 and variance 1/a and so we may, 
for large a, view it as a perturbation of yU°°. Indeed, by using our explicit formula 
for A^a together with the large-order expansions in Abramowitz & Stcgun (1964), 
§9.7, we find 

1 l + 8t (\ 
a2(l + 4i)^ \a 

Likewise, setting 



Apa (t) ^ A^oo (t) + - + O — ) asa^oo, a e Z+ . 



/3 := arccos , < A < 4 , 

and using the large-order expansions in Abramowitz & Stegun (1964), §9.3, we 
obtain, for a G Z+, 

q{\) ^ (A) - 1-^^^ (13 + 38 cot^ /? + 25 cot^ + O ( 
327rsin (3 ' 

This expansion breaks down at A = 4; as a increases, g{X) diverges to infinity there 
but tends to zero exponentially fast for A > 4. 

The following computations were performed in multiple-precison floating-point 
arithmetic with the MAPLE software package; the eigenvalues and eigenvectors of the 
matrix ^„ were calculated by using the Eigenvals function, which implements 
the QR algorithm. Figure [1] corresponds to the case where = with a = 8; it 
illustrates the convergence of the counting measure Nn to the integrated density of 
states N as n increases — thus confirming the validity of our formula for N. 
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(i) (ii) 

Figure 1. The counting measure Nn (solid line) for a particular 
realisation of the sequence s when ^ = with a = 8: (i) n = 128 
and (ii) n — 256. For comparison, points corresponding to values 
of the integrated density of states N are also shown. 



While the Lyapunov exponent and the density of states are non-random, the 
measure a is random. We note that a°° is absolutely continuous whereas, for every 
a e Z+, almost every realisation of a is singular. Figure [2] shows the approximation 

/ ^T„(dA'), 
Jo 

of the integrated measure for particular realisations corresponding to a = 8 and 
a = 64. Here (t„ is the discrete measure defined by the quadrature formula (jl.lOp . 
and n = 256. 
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Appendix A. Proof of proposition 11.11 

As mentioned in the introduction, the theory of products of random matrices is 
a convenient means of deducing the uniqueness of the invariant measure, as well as 
the positivity of the real part of the complex Lyapunov exponent. For this purpose, 
we shall have to deal with products of 2 x 2 matrices with real or complex entries; 
so, in the following, K will stand for cither K or C. Set 

K KU {oo} . 

We define an equivalence relation in the set of nonzero vectors in via 

") ^ ("') ^ ^ ^"'^ (") ^ 

The set of the equivalence classes is called the projective space P(]K^). Let 

e P (K^) . 



We shall identify this equivalence class with 

\ u/v ifuT^O 



V 



oo otherwise 



The results of Furstenberg & Kesten (1960) and Furstenberg (1963) concern the 
typical asymptotic behaviour of the product of independent, identically distributed 
random elements of some group acting on a compact topological space. In our 
particular context, the relevant group is the subgroup of GL (2,K^) consisting of 
2x2 matrices with determinant ±1, and the topological space is P(K^). The 
invertible matrices 



(A.l) A = 



a b 
c d 



are drawn at random from a distribution which we shall denote by jl. The action 
of the matrix A on the projective space can be expressed as 

1 



(A.2) ^.1 



where T : K ^ K is the linear fractional transformation defined by 

2^ if z e K and cz + d 7^ 
OO if c 7^ and z — —d/c 



(A.3) T{z)^V[A 1 , , 

\ J J o./c it z = oo and c (J 

oo if c — and z = oo 
Thus, we have an obvious connection between products 

lAn '■— An ■ ■ ■ A2A1 
of the An and the Markov chain Z such that 

(A. 4) Zo = z, Z„+i = Tn [z^ . 
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Indeed, 
(A.5) 



We are interested in the rate of growth of this product as n — s- cx); this is quantified 
by the number 



(A.6) 



7^ hm -E(ln|ZY„|), 

n— >oo Tl 



where | • | denotes some matrix norm. This hmit exists whenever 
(A.7) E(log+ |yl|) < cx). 

In the htcrature on products of random matrices, the name "Lyapunov exponent" 
is usuaUy reserved for 7^ ; as we shall see, it is in fact the real part of the complex 
Lyapunov exponent introduced earlier. 

The following specialisation of Furstenberg's theorem will be the most useful for 
our purpose: 

Theorem 4. Let the An be independent 2x2 random matrices with determinant 
±1 and a common distribution jl. Suppose that there is no measure on P (iK^) that 
is invariant under the action of the smallest subgroup o/GL(2,K^) generated by the 
support of jl. Then 

(1) there exists a unique probability measure Vjx on P (K^) that is invariant 
under the action of matrices in the support of fi. 

(2) Let z G K. Then, for almost every realisation of the sequence {A^l\n&i, 



lim — In 

n^ca n 



7/1- 



(3) // condition IjA.!^ holds, then the number 7^ is strictly positive and is given 
by the formula 



7/i 



In 



A 



GL(2,K) 



fL (dA) Vf, d 



Proof. See Bougerol & Lacroix (1985), Part A, theorem 4.4 on p. 32. The last 
formula appears in Part A, theorem 3.6 on p. 27 of the same reference. □ 

Corollary A.l. Let t E C\IR_. Let the s„ be independent random variables in 
with a common distribution /j, whose support contains at least two points. Set 



1 



Then fi induces a probability distribution jl on GL(2,C) such that the hypothesis of 
theorem^is fulfilled, and the image of the invariant measure vp, under the map 
V is the unique stationary distribution for the Markov chain 7i. Furthermore, if we 
also assume that 



I s'^ A^(ds) < 00 
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for some e > 0, then we have the formula 

7a = - / ^n\z\i^^{dz). 



Proof. Same as that of Carmona & Lacroix (1990), corollary IV. 4. 26. □ 

Corollary A. 2. Let u„ be the sequence defined by the recurrence relation I11.16\) 
with uq 7^ 0. Under the same hypothesis on ii, for every t e C\M_, for almost every 
realisation of the sequence s, we have 

lim 1^=A^(0 

n — *oo Tl 

independently of the starting value uq. 
Proof. This follows from the equality 

In Un In Uq 



1 " 



n n n 

i=i 

the definition of A^, and the uniqueness of the stationary measure for the Markov 
chain Z. □ 

For reasons that will become clear when we discuss the nature of the spectrum 
of the measure a, we also need to consider the case where 

t^-x + iO±, x>0. 

In this case, 

'0 1 



1 



and so the recurrence defining the Markov chain Z is 

1 



Zn+l — 



Once again, theorem 3] implies the existence of a unique stationary measure fi^. We 
remark that this stationary measure is concentrated on the imaginary axis. Indeed, 
define ¥„ via 



Then 

Y - 



The resulting Markov chain Y — {Yq ,Yi , . . .) corresponds to the product of the 
real (Schrodinger) matrices 

By applying theorem [3] with K = R, we deduce that Y has a unique stationary 
measure supported on R. The fact that the measure p^i, must be concentrated on 
the imaginary axis then follows by uniqueness. 
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Corollary A. 3. Let Un be the sequence defined by the recurrence relation il.l6\) 
with uo ^ and t = — x + iO±. Under the same hypothesis on fi, for every x € R+, 
for almost every realisation of the sequence s, we have 

lim i^i!^ = A J-x + iO±) 

n — *OD Tl 

independently of the starting value Uq. 

Finally, the positivity of 7^, asserted in theorem SI leads to the 
Corollary A. 4. Under the same hypothesis on 11, for every t G C\M-, 

Re[A^(0] > 0. 

Also, for every x G M_|_, 

Re[A^(-x + iO±)] > 0. 



Proof. Take the real part in equation (|1.20p and use corollarv lA.il . □ 

The statements made in corollaries IA.mA.4l are of the form: 

Let t be fixed, then for almost every realisation of the sequence s, 
etc. 

But the proof of proposition 11.11 follows easily from them by a well-known argument 
based on the use of Fubini's theorem; see, for instance, Ishii (1973). 

Appendix B. The complex Lyapunov exponent for the gamma 

distribution 

In this appendix, we derive the formula (jl.26p for the complex Lyapunov expo- 
nent when /i is the gamma distribution. For this purpose, it is convenient to adopt 
the notation used in Marklof et al. (2005); so we set 

aigt 



and 



e [-7r/2, 7r/2] . 



Then 



aie" 



1 



a2e" 



ase" 



The random variable Z takes values in the set 

S'a := {z G C : |argz| < \a\} . 

In Marklof et al. (2005), we showed that, if the a„ are gamma-distributed with 
parameters p and s, then the probability density function of Z is given explicitly 

by 



(B.l) fo.{z) 



sin(2|Q;|) 



1 



|2Xp(2e'"/s)r r2 sm\a + e) 



sm a 



sin(Q; -|- I 



X exp 



sin(2Q;) 



r sin(Q! 



sin(Q; -I- 9) 
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where r = |z| and 9 = aigz. This is vaUd for \a\ < 7r/2 which is all wo need 
for the purpose of calculating the complex Lyapunov exponent, since the cases 
a = ±7r/2 may be obtained by letting a tend to the appropriate limit. In Marklof 
et al. (2005), we derived the formula 



-I 



ln|^;|/Q,(^) dz = Re 



So there only remains to show that 

(B.2) - / In (e^'^^'s^) foc{z) dz = ilm 



9pKp (f e--] 
K, (f e-i«) 



We shall need the 

Lemma B.l. Let iv and W be two complex numbers with positive real part. Then 

2v 



In- 



Kp{w)dpKp{W)=-l- [ e-2p« / fu , 

+ 2wW cosh(2u) + W\ dv 



X exp 



4v 



— du . 

V 



Proof. The Bessel function of the second kind has the integral representation (see 
Watson 1966, §6.22) 



1 f 

Kp{z) = -J e-^^°^^^-f^da;. 



Hence 



Kp{W)dpKp{w) = / / e:>ip{— W cosh y — py — w cosh X — px} xdxdy . 



We make the change of variables 

x = X + Y and y = X-Y . 
Then, after some re- arrangement, 



Kp{W)dpKp{w)=~ I e 



-2pX 



X / exp < -e 6 + YfdYdX . 



Make the substitution 

in the second integral. Then 

Kp{W) dpKpiw) = I e 



2 2 



V = e 



-2pX 



r / 

< / exp l—v — 



+ 2wWcosh{2X)+w^ 



X + \n: 



2v 



dv 



dX. 



The desired result follows after we set X = u and take the imaginary part. □ 
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Returning to the proof of equation (jB.2|) . let us write 

1 



(B.3) 
where 
X := Im 



Im 



In 



In(e-'^'-g^) Uiz)dz 



sin(2|a|) 



4|X,(|e--) 



|2 ^ ' 



a Jo 



sin (a + i 



sin(Q; + I 



X exp 

Replace the variable r by 
and then the variable 9 by 



sin(2a) 



1 



r sm a 



sin(Q; 



drde. 



t = 



r sin(2a) 
s sin(a + 

sin(a — 6) 



sin(Q! + I 



A straightforward calculation leads to 
(B.4) 



I = Im 



In 



tP- 



exp < —V 



s 4u 



where 



Set 



(p{t) = t + - + 2 cos a . 



t = e 



-2u 



Then equation (IB.4[) becomes 



(B.5) X = -2 Im / In (e"+'" + e"""'") e^^P" 



exp < —V 



4 e 



-i2Q 



2 cosh(2M) + e" 



i2a 



dv 

— du . 

V 



We use lemma IBTT] with = w = 2c '"/s to obtain 



(B.6) X=-4Im 



K, [ ^c- 



dpKp I -e- 



4Im 



The result is then an immediate consequence of equation (|B.3[) . 

Appendix C. A result of Goldsheid & Khoruzhenko (2005) 

In order to deduce the existence of the integrated density of states from proposi- 
tion l2.1i we shall need some technical results contained in Goldsheid & Khoruzhenko 
(2005). 

Let {^n}ngz+ be a deterministic sequence of square matrices of increasing di- 
mension n, and let 

Pn{z) = - In |det {x/n - zJn)\ = I \n\w - z\ Kn{dw) , 
n Jc 
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where J'n is the n x n identity matrix and 

is the normalised eigenvalue counting measure of Define 
(C.l) Tfl:=Ttar/ ln|w|K„(dw), R>1. 

Proposition C.l. Assume that there is a function p : C [—00, 00) such that 
Pn{z) > p{z) for Lehesgue-almost every z G C . 

n — >oo 

If Ti < 00, then 

(1) p is locally integrable. 

(2) The measure 

1 ^ 
K := — An 
2tt ^ 

is a probability measure. 

(3) We have 



/ In \w\ K{dw) < Ti < 00 

J\w\>l 



l\w\>l 

and the sequence {Kn}nez^ converges weakly to n. 
Proof. See Goldsheid & Khoruzhcnko (2005), proposition 1.3. □ 



To apply this result in our context, we set £/n = ^n- Proposition [2A] then takes 
care of the first assumption. There remains to show the finiteness of t\. Goldsheid 
& Khoruzhenko remark that the following inequalities hold: 

Ti < TtaT-^trln (A + A A*) 

and 

^ n— 1 

(C.2) _trln(A + AA*)<^El^(l + /3|^^l') ' 

for some positive constants a and (i independent of n and r^, the jth row of 
Now, 



^hl_^+vl^hl 



1 1/1 1 



+ 



52n+2'52n+lS2n+3 sL+l V S2n+2 S2ri+3 / si„'S2n-lS2n+l 

and, so by repeated use of the elementary inequality 



we obtain readily 



.n 5 
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It follows that 



/ 25 

In (1 + /?|r„|2) < In (^1 + _/3niax {s2;f_i, . . . , 



< In f 1 + ^l3s^,Uj + • • • + In (^1 + f /3s2-„\3 

We deduce from the ergodic theorem that the right-hand side of equation (jC.2[) is 
bounded independently of n provided that 

This last inequality follows easily from our hypothesis that 

/ |ln s| /i(ds) < CX3 and / s^^(ds)<oo. 
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